Solid-solid volume collapse transitions are zeroth order 
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We present an exactly solvable non-linear elastic model of a volume collapse transition in an 
isotropic solid. Integrity of the lattice through the transition leads to an infinite-range density- 
density interaction, which drives classical critical behavior. Nucleation is forbidden within a pres- 
sure window leading to intrinsic hysteresis and an unavoidable discontinuity of the thermodynamic 
potential (zeroth order transition). The window shrinks with increasing temperature ending at a 
critical point at a temperature related to the shear modulus. Mixed phases behave non-extensively 
and show negative compressibility. We discuss the implications for Ce, SmS, and related systems. 

PACS numbers: 64.70.Kb, 68.35.Rh, 71.28.+d, 71.10.Hf 



Room temperature 7 Cerium shows a well known tran- 
sition to isostructural a Cerium at ^ O.SGPa with a 
17% reduction of the volume. Similar (though not neces- 
sarily isostructural) transitions occur in other /-electron 
systems 0,0- Volume-collapse transitions were also pre- 
dicted to occur in colloidal systems Q and close to first 
order electronic transitions when mesoscopic inhomoge- 
ncous electronic states are frustrated by the long-range 
Coulomb interaction^. 

Consider one of these systems which is driven across 
the volume-collapse transition by adiabatically changing 
the average volume. In analogy with the van der Waals 
liquid-gas transition one expects a range of average vol- 
ume values for which there is phase coexistence at a con- 
stant equilibrium pressure. This presents an interesting 
nucleation problem. Within purely energetic considera- 
tions the mixed phase consists of a region of the sample 
with the small unit cell volume and a region with the 
large unit cell volume separated by a flat surface so as 
to minimize the interface boundary energy. The lattice 
mismatch at the interface can be seen as a line of dislo- 
cations. However dislocations are very difficult to nucle- 
ate in a solid since they require the breaking of chemical 
bonds. Indeed it is customary to assume on solid-solid 
transitions that the lattice integrity is preserved 0, 0,0] ■ 
Under this hypothesis the thermodynamic ground state 
is frustrated from the onset, and the question arises of 
how the system minimizes its free energy, preserving the 
integrity constraints. To answer this question we solve a 
non-linear elastic model of a volume-collapse transition 
in an isotropic solid assuming lattice integrity (no dislo- 
cations). The latter condition leads to the well known 
St. Venant compatibility constraints among the differ- 
ent components of the strain tensor and allows the free 
energy to be expressed as a function of the order parame- 
ter (OP) alone, with long-range interactions 0,0,0]. This 
kind of approach has been shown to lead to anomalous 
nucleation in volume conserving transitions'^. 



We find that compatibility constraints induce a 
distance-independent density-density interaction and 
thus mean-field (MF) theory becomes exact and the 
model results to be trivially solvable at all temperatures 
in the absence of short-range interactions. Critical ex- 
ponents are MF like even in the presence of short-range 
interactions. In the mixed phase the system behaves non 
extensively and violates usual stability criteria for ex- 
tensive systems. In particular, it possesses negative com- 
pressibility. For an experiment in which the pressure acts 
as the control variable we show that usual nucleation is 
forbidden within a finite pressure window. This pressure 
window implies an intrinsic hysteresis loop in a pressure 
controlled experiment. Outside the pressure window the 
transition can occur, but only irreversibly with an in- 
trinsic jump of the thermodynamic potential. Since the 
discontinuity occurs in the thermodynamic potential it- 
self rather than in its first derivative the transition can be 
called of zeroth order (instead of first) as done with sim- 
ilar phenomena in gravitational systems 0. The extent 
of the hysteresis loop reduces as a function of tempera- 
ture, vanishing at a critical point. The temperature of 
the critical point is related to the shear modulus of the 
material. 

Model The OP is taken to be the dilation strain 
ei = Y^'^^i £u / Vd where ejj are the components of the 
infinitesimal strain tensor [Sj e and d is the dimensionality 
of the system. Our starting point is a coarse grained lO^ 
elastic Hamiltonian H = Hi + Hh with Hi a non-linear 
Ginzgurg-Landau like double-well term: 



Hi^ j d'^x (^aei - ^-el + '-e{ + f iVeij^) 



(1) 



The remaining d{d -t- l)/2 — 1 independent components 
of the strain tensor give the harmonic contribution: 
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FIG. 1: (a) (^-dependent common tangent construction. The 
cases (j) = —0.5 and 0.2 are illustrated (m = 1). The tri- 
angles indicate the minimum value of h attainable for the 
corresponding values of (j). (b) Analytical results for h{(j}) as a 
function of (j) for 7 = and different values of m, as indicated. 
For m > 2 the solution is h{(j>) — ho{(j}). 



where fi is the shear modulus of the material, and the 
(traceless) 'accommodation strains' tensor = — 
SijCi/Vd has been defined. 

Zero temperature theory In the absence of disloca- 
tions and in order to be physical the strain tensor must 
fulfill the St. Venant compatibility conditionj^]: 



curl curl e = 



(3) 



i/, = (d-l)/i-L^ei(k)ei(-k). 



Here L is the linear dimension of the system. A crucial 
point is that the k = term is excluded from the sum. 
We can extend the sum to all k by explicitly subtracting 
the missing term. We obtain in real space: 



H,J[{d-l)^I] - / d''x[e,{^)-e 



(4) 



with ei = / d^'xei/L'*. The accommodation strains in- 
duce a local term that increases the stiffness at all wave 
vectors, and a distance independent "ferro" like interac- 
tion in ei that decreases the uniform stiffness alone, i.e., 
the stiffness at finite wave vector is larger than at zero 
wave vector, a well known result that holds also for linear 
elastic solids [foj. A similar result has been obtained in 
d = 2 in the context of the dislocation mediated melting 
problem close to an isostructural critical point 0]. 

We define (/) = ei^Jc/b, m = 4(d— and 7 = 2k/6. 

Up to an additive constant the coarse grained energy 
density h [from Eqs. and in units of the height of 
the double well barrier [&^/(4c)] reads: 



(5) 



which leads to constraints between ei and e^j . A uniform [1 
strain satisfies Eq. trivially and therefore it is uncon- 
strained. Thus the k = term in the Fourier decompo- 
sition of iJ/i, corresponding to uniform accommodation 
strains, can be minimized by choosing eij(k = 0) = 0. 
For k 7^ the accommodation strains can be elimi- 
nated in favor of ei by minimizing Hf^ with respect to 
Zij , and enforcing the constraints in Eq. (|3Jl by Lagrange 
multipliers!^. The resulting equilibrium equations to- 
gether with Eq. ^ allow Eq. (01 to be written as a func- 
tion of ei alone: 



ei(x)ei(x') 



with /io(0) = (1 — 0^)^ being the double well in the di- 
mensionless volume 0. The term linear in the OP has 
been eliminated by redefining the origin of pressure. The 
parameter 7 fixes the width of the interface between the 
low and high volume phases in a mixed state. 

We first discuss the 7 = case and generalize for a 
finite 7 below. The energy density Eq. © couples the 
values of (j) at different spatial positions only through 
0, namely it corresponds effectively to a MF. We stress 
however, that this is the exact energy of the model. 

Mixed phase solutions for a constrained total vol- 
ume obey the familiar common tangent construction: 
dh{(j),4))/d(l)\4>=4,_ = dh{(l),4>)/d(j}\^= ^^ = [h{(j )^,$) - 
(f>)]/{(j)- — (/)+) where (j)± = ± •^/(T^— to72)' and the 
+ (— ) corresponds to the OP of the expanded (collapsed) 
phase. The only difference with a conventional phase 
separation computation is that the construction is (j)- 
dependent because h depends on the global variable (j) 
see Fig. nja)]. This dependence makes the total energy 
of the mixture non additive. In Fig. ^b) we report the 
resulting average energy density h{(f)). For m > 2 the 
mixed solution does not exist and the system remains al- 
ways uniform. We have checked these results with direct 
analytical computations assuming specific shapes for the 
coexisting phases and different boundary conditions. In 
all cases the maximum shear strain is of the same order 
as the relative volume mismatch. As long as the volume 
mismatch is small (a few percent) the shear strain will 
behave elastically. For large volume mismatch glide of 
atomic planes may take place and our theory becomes ap- 
proximate. In any case the transition will be irreversible. 

For a pressure controlled experiment the physical pres- 
sure IS given by P = p^/WJ{16cd) - a/Vd with p the 
dimensionless pressure: p = —dh{(f>)/d(f> at T = 0. In the 
mixed phase the free energy has negative curvature so 
the compressibility is negative [Fig. Wi^)] ■ This thermo- 
dynamic state is stable at constant volume with the pined 
boundary conditions of Ref. • This violation of usual 
stability criteria is common in non-extensive systems |^. 

In Fig. 121 we show the values of cj) and the Gibbs free 
energy density g{p) = /i[0(p)] + p(j){p) as a function of p 
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for two values of m. Controlling the total volume the sys- 
tem would follow the continuous line. The straight por- 
tion of this solution in Fig. |2Ia) corresponds to the mixed 
phase. Consider now a pressure controlled experiment for 
m < 4/3 [Fig. [2Ia) and (c)]. Starting from the high vol- 
ume phase at low pressure and increasing p the system 
remains uniform until p^, given by the point in which 
the mixed solution line intersects the uniform solution. 
At this point the small volume phase nucleates sponta- 
neously since 7 = 0. It is easy to show that the mixed 
state is mechanically unstable at constant pressure due 
to the negative compressibility, so the volume collapses 
and the system jumps to the lower branch of the curve 
following the vertical line. Decreasing the pressure the 
process is reverted and the volume suddenly expands at 
p^{= —p^) producing the hysteresis loop indicated in the 
figure. A conventional first-order phase transition driven 
by pressure at a very low rate (so that nucleation occurs) 
would follow a vertical line at p = in the top panel. 
This is not possible here because the infinite range inter- 
action makes the critical drop radius diverge, and the 
nucleation energy-barrier to scale as the volume of the 
system. 

In the case m > 4/3 the compressibility of the uni- 
form solution diverges at the spinodal points [full dots in 
Fig. 111(b) and (d)]. The transition proceeds by a k = 
instability at p^ — —p^ = 8/(3V3), independent of m 
and with no energy barrier. The lower panels in Fig. 12 
show that the transition involves an unavoidable jump 
in the thermodynamic potential for all m ^ justifying 
calling this transition 'zeroth order' 0. 

Finite temperature theory Using MF theory it is 
straightforward to obtain the exact partition function for 
7 = 0, and the exact (f> ~ P dependence at different tem- 
peratures. In Fig. 121 we show the isotherms in the (j) ^ p 
plane for m = 0.6 for different values of a dimensionless 
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FIG. 2: Mean volume (j> and Gibbs free energy density for 
m — 0.6 (left) and m > 2 (right). Dotted lines are the uni- 
form solution. Continuous lines are the solutions that min- 
imize Eq. (jSJ for each value of 4>. The arrows and dashed 
lines indicate the hysteresis cycle in pressure controlled ex- 
periments. Full dots correspond to p\ p^. 



temperature t (units are restored below). The pressure 
width (Ap) and volume width (A^) of the hysteresis loop 
is reduced at finite temperatures up to the point where it 
vanishes, which determines the critical point coordinates 
{pc = Ojtcim)). The pressure width as a function of 
temperature determines a wedge of forbidden nucleation 
[shown in Fig. j^b) for different values of m] around the 
critical pressure pc — 0. 

For m < 4/3 the singularities dX p — of the 

T = isotherm become rounded at finite temperature 
[Fig. Oa)] due to the thermally induced spatial fluctu- 
ations of the OP. This implies that the compressibility 
will show precursor effects for all m > 0, diverging at the 
transition. 

The solution in the case of finite 7 and m is given im- 
plicitly by = 00 (p + ™0) where 4>q {p) is the equation 
of state of the same model except that m = 0, i.e., the 
well studied (f>^-m.ode\ which has the universality class of 
the short-range Ising model|0. If m > 4/3 the quali- 
tative description of the forbidden nucleation region are 
completely equivalent to the 7 = 6 case. For to < 4/3 a 
region of normal nucleation and growth appears around 
the wedge of forbidden nucleation at low temperatures 
and the divergence of the compressibility gets a cutoff. 
These effects never reach the critical point. Both for 
7 = and 7 > we find that critical exponents are given 
exactly by MF. This is also in agreement with the d — 2 
theoretical analysis of Ref. ■ 

Comparison with experiments A wedge where nucle- 
ation does not occur is observed in Celll Q, SmS0 
and related systems. Interestingly it is also observed in 
a volume-collapse transition in amorphous ice|l4|. How- 
ever without detailed studies it is difficult to be sure that 
it does not arise from usual hysteresis due to a finite 
driving rate, as occurring in other first order transitions. 
Other indicators however, point more strongly to the ap- 
plicability of our theory. 

In agreement with the present theory MF critical ex- 
ponents have been observed^. The connection between 
MF critical behavior and long-range strains has long been 
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FIG. 3: (a) (j> vs. p for m = 0.6 and different values of the 
dimensionless temperature t, as indicated, (b) t — p phase 
diagram for 7 = and various m values (indicated in the 
figure). The interior of the wedge (for each value of m) is a 
region in which nucleation is forbidden. 
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suspected 121 ■ 

In contrast with usual first-order phase transitions, 
where hysteresis is not intrinsic, one can define here a 
critical exponent for the pressure width of the hysteresis 
loop: Ap oc |i — icP^^- We are not aware of detailed mea- 
sures of this critical exponent, however the curvature of 
the boundaries found in Ref. flJl for Gd doped SmS are 
compatible with this exponent. We hope our work will 
stimulate further experimental work. 

Also in contrast with usual first-order transitions, criti- 
cal behavior (although classic like) can be easily observed 
far from the critical point. The divergence of the com- 
pressibility discussed above implies that the bulk modu- 
lus vanishes as Ip — and jumps discontinuously to a 
finite value for p > . Decreasing pressure the jump and 
the square root singularity change side and occur at p'^ . 
Similar behavior is found as a function of temperature. 

An anomaly measured in YbInCu4 on cooling resem- 
bles our prediction when plotted in an appropriate scale 
(see inset of Fig. 4 in Ref. ^3)- This also explains the 
appearance of precursor effects which have been very puz- 
zling in the mixed valence literature and have led some 
authors to interpret the transition as second order [T^. 

In mixed valence systems the double well is attributed 
to the anomalous contribution of /-electrons to the crys- 
tal binding [l^. For a realistic description it may be nec- 
essary to take the parameters in Eq. Q to depend on 
temperature to take into account temperature effects of 
the /-electrons. In fact a trivial linear and quadratic tem- 
perature dependence of a fits our symmetric phase dia- 
gram to the asymmetric forms of CeP,Q and SmSf^Hsj. 

The temperature dependence of a does not affect the 
critical temperature, but b does. Indeed in current ex- 
planations of the critical point in /-electron systems this 
electronic effect is crucial whereas the shear rigidity does 
not play a relevant role|0. To check the relevance of 
the temperature scale set by the shear modulus on de- 
termining real values of Tc we take the opposite point of 
view, i.e., we neglect the temperature dependence of the 
coefficients in Eq. ^ and take into account the effect 
of the shear rigidity alone. Restoring dimensions, the 
critical temperature is given by: ksTc = tc(jn)vob'^ / (Ac) 
where is the minimum volume in the system allowed 
to fluctuate independently. For the case 7 = this vol- 
ume should be assimilated to the atomic volume. In the 
case 7 > 0, this volume can be greater depending on the 
strength of the short distance interactions parametrized 
by 7. In any case we expect it to be of the order of a 
few atomic volumes. The parameters 6, c and m can 
be fixed by fitting the observed pressure and volume 
with of the hysteresis loop at a given temperature and 
from the measured value of the shear modulus P, ll^ . 
Within the experimental uncertainties this procedure 
gives the right value of for Ce (- 550k[iI Q]) and 
SmS (~ 750X0, 01) if we assume vq ~ 2vf where Vf 
is the volume associated with one formula unit. Since 



the relation between vq and Vf go beyond our model we 
cannot exclude temperature effects on the coefficients of 
Eq. (0) contributing to however the fact that we 

obtain the right order of magnitude with a vq on the or- 
der of Vf implies that the scale set by the shear modulus 
can not be neglected in the determination of Tc- 

In conclusion we have presented an exactly solvable 
model of an isostructural volume collapse transition, the 
solid-state analog of the van der Waals transition. Unlike 
the fluid counterpart the thermodynamics is very anoma- 
lous due to the presence of inflnite range interactions. In 
particular, the transition is intrinsically irreversible with 
a discontinuity of the thermodynamic potential making 
of it a laboratory example of a class of phase transitions, 
namely zero order transitions, flrst discussed in the con- 
text of gravitational systems |9j. 

Despite the simplicity of our model, which assumes an 
isotropic solid, we believe the present features are quite 
robust. Indeed our flndings are in good agreement with 
the behavior of real anisotropic svstemsflTO. Iisl ITsf . In 
these mixed valence systems the unexpected thermody- 
namics is a consequence of the interplay between st rong 
correlations, determining the double well potential[la|, 
and long-range strains. The present approach can be 
seen as a flrst step to understand this interplay in more 
complex non-isostructural transitions like the magnetore- 
sistant manganites'i'l which also show volume mismatch 
among different phases and irreversibility 18J. 
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